Investigating the importance of left atrial compliance on fluid dynamics in a novel mock circulatory loop.

The left atrium (LA) hemodynamic indices hold prognostic value in various cardiac diseases and disorders. To understand the mechanisms of these conditions and to assess the performance of cardiac devices and interventions, in vitro models can be used to replicate the complex physiological interplay between the pulmonary veins, LA, and left ventricle. In this study, a comprehensive and adaptable in vitro model was created. The model includes a flexible LA made from silicone and allows distinct control over the systolic and diastolic functions of both the LA and left ventricle. The LA was mechanically matched with porcine LAs through expansion tests. Fluid dynamic measures were validated against the literature and pulmonary venous flows recorded on five healthy individuals using magnetic resonance flow imaging. Furthermore, the fluid dynamic measures were also used to construct LA pressure-volume loops. The in vitro pressure and flow recordings expressed a high resemblance to physiological waveforms. By decreasing the compliance of the LA, the model behaved realistically, elevating the a- and v-wave peaks of the LA pressure from 12 to 19 mmHg and 22 to 26 mmHg, respectively, while reducing the S/D ratio of the pulmonary venous flowrate from 1.5 to 0.3. This model provides a realistic platform and framework for developing and evaluating left heart procedures and interventions.

Novel Robust Needle Tip Design Enables Needle Reuse and Reduced Skin Trauma in Combination With Autoinjector Needle Shields.

BACKGROUND: Pen needles and autoinjectors are necessary for millions of patients needing injectable drug treatment but pose economic and environmental burdens. A durable device with a multiuse needle could reduce cost and improve user experience. This study explores a novel robust needle tip (EXP) designed for multiple uses and durability against hooking. METHOD: Needle robustness was investigated through a structural analysis. Furthermore, EXP and control needles (NF30, NF28) were evaluated in an in-vivo porcine model as pen needles or embedded in autoinjectors to study the resulting increase in skin blood perfusion (SBP). The SBP was assessed by laser speckle contrast analysis (LASCA) of 192 randomized and blinded needle insertions. RESULTS: Forming a 33 µm hook against a hard surface requires 0.92 N for the NF30 control needle and 5.38 N for EXP. The EXP did not induce more tissue trauma than the NF30. There was a positive relation between needle diameter and SBP (P < .05). Furthermore, the presence of an autoinjector shield and applied force of 10 N was found to significantly reduce SBP for worn EXP needles (P < .05) compared to insertions without autoinjector shield. CONCLUSIONS: The investigated robust needle EXP is on par with the single-use needle NF30 in terms of tissue trauma, which is further reduced by combining the needle with a needle shield. These results should encourage the innovation and development of durable, reusable injection systems with pharmacoeconomic and environmental value and a simplified and enhanced user experience for patients.

Intra-Cardiac Flow from Geometry Prescribed Computational Fluid Dynamics: Comparison with Ultrasound Vector Flow Imaging.

PURPOSE: This paper investigates the accuracy of blood flow velocities simulated from a geometry prescribed computational fluid dynamics (CFD) pipeline by applying it to a dynamic heart phantom. The CFD flow patterns are compared to a direct flow measurement by ultrasound vector flow imaging (VFI). The hypothesis is that the simulated velocity magnitudes are within one standard deviation of the measured velocities. METHODS: The CFD pipeline uses computed tomography angiography (CTA) images with 20 volumes per cardiac cycle as geometry input. Fluid domain movement is prescribed from volumetric image registration using the CTA image data. Inlet and outlet conditions are defined by the experimental setup. VFI is systematically measured in parallel planes, and compared to the corresponding planes in the simulated time dependent three dimensional fluid velocity field. RESULTS: The measured VFI and simulated CFD have similar flow patterns when compared qualitatively. A quantitative comparison of the velocity magnitude is also performed at specific regions of interest. These are evaluated at 11 non-overlapping time bins and compared by linear regression giving R2 = 0.809, SD = 0.060 m/s, intercept = - 0.039 m/s, and slope = 1.09. Excluding an outlier at the inlet, the correspondence between CFD and VFI improves to: R2 = 0.823, SD = 0.048 m/s, intercept = -0.030 m/s, and slope = 1.01. CONCLUSION: The direct comparison of flow patterns shows that the proposed CFD pipeline provide realistic flow patterns in a well-controlled experimental setup. The demanded accuracy is obtained close to the inlet and outlet, but not in locations far from these.

Precise Estimation of Intravascular Pressure Gradients.

This study presents a method for noninvasive pressure gradient estimation, which allows the detection of small pressure differences with higher precision compared to invasive catheters. It combines a new method for estimating the temporal acceleration of the flowing blood with the Navier-Stokes equation. The acceleration estimation is based on a double cross-correlation approach, which is hypothesized to minimize the influence of noise. Data are acquired using a 256-element, 6.5-MHz GE L3-12-D linear array transducer connected to a Verasonics research scanner. A synthetic aperture (SA) interleaved sequence with 2 ×12 virtual sources evenly distributed over the aperture and permuted in emission order is used in combination with recursive imaging. This enables a temporal resolution between correlation frames equal to the pulse repetition time at a frame rate of half the pulse repetition frequency. The accuracy of the method is evaluated against a computational fluid dynamic simulation. Here, the estimated total pressure difference complies with the CFD reference pressure difference, which yields an R -square of 0.985 and an RMSE of 3.03 Pa. The precision of the method is tested on experimental data, measured on a carotid phantom of the common carotid artery. The volume profile used during measurement was set to mimic flow in the carotid artery with a peak flow rate of 12.9 mL/s. The experimental setup showed that the measured pressure difference changes from -59.4 to 31 Pa throughout a single pulse cycle. This was estimated with a precision of 5.44% (3.22 Pa) across ten pulse cycles. The method was also compared to invasive catheter measurements in a phantom with a 60% cross-sectional area reduction. The ultrasound method detected a maximum pressure difference of 72.3 Pa with a precision of 3.3% (2.22 Pa). The catheters measured a maximum pressure difference of 105 Pa with a precision of 11.2% (11.4 Pa). This was measured over the same constriction and with a peak flow rate of 12.9 mL/s. The double cross-correlation approach revealed no improvement compared to a normal differential operator. The method's strength, thus, lies primarily in the ultrasound sequence, which allows precise and accurate velocity estimations, at which acceleration and pressure differences can be acquired.

Anatomic and Functional Imaging Using Row-Column Arrays.

Row-column (RC) arrays have the potential to yield full 3-D ultrasound imaging with a greatly reduced number of elements compared to fully populated arrays. They, however, have several challenges due to their special geometry. This review article summarizes the current literature for RC imaging and demonstrates that full anatomic and functional imaging can attain a high quality using synthetic aperture (SA) sequences and modified delay-and-sum beamforming. Resolution can approach the diffraction limit with an isotropic resolution of half a wavelength with low sidelobe levels, and the field of view can be expanded by using convex or lensed RC probes. GPU beamforming allows for three orthogonal planes to be beamformed at 30 Hz, providing near real-time imaging ideal for positioning the probe and improving the operator's workflow. Functional imaging is also attainable using transverse oscillation and dedicated SA sequence for tensor velocity imaging for revealing the full 3-D velocity vector as a function of spatial position and time for both blood velocity and tissue motion estimation. Using RC arrays with commercial contrast agents can reveal super-resolution imaging (SRI) with isotropic resolution below [Formula: see text]. RC arrays can, thus, yield full 3-D imaging at high resolution, contrast, and volumetric rates for both anatomic and functional imaging with the same number of receive channels as current commercial 1-D arrays.

Performance Assessment of Row-Column Transverse Oscillation Tensor Velocity Imaging Using Computational Fluid Dynamics Simulation of Carotid Bifurcation Flow.

In this work, the accuracy of row-column tensor velocity imaging (TVI), i.e., 3-D vector flow imaging (VFI) in 3-D space over time, is quantified on a complex, clinically relevant flow. The quantification is achieved by transferring the flow simulated using computational fluid dynamics (CFD) to a Field II simulation environment, and this allows for a direct comparison between the actual and estimated velocities. The carotid bifurcation flow simulations were performed with a peak inlet velocity of 80 cm/s, nonrigid vessel walls, and a flow cycle duration of 1.2 s. The flow was simulated from two observation angles, and it was acquired using a 3-MHz 62+62 row-column addressed array (RCA) at a pulse repetition frequency ( fprf ) of 10 and 20 kHz. The tensor velocities were obtained at a frame rate of 208.3 Hz, at fprf = 10 kHz , and the results from two velocity estimators were compared. The two estimators were the directional transverse oscillation (TO) cross correlation estimator and the proposed autocorrelation estimator. Linear regression between the actual and estimated velocity components yielded, for the cross correlation estimator, an R 2 value in the range of 0.89-0.91, 0.46-0.77, and 0.91-0.97 for the x -, y -, and z -components, and 0.87-0.89, 0.40-0.83, and 0.91-0.96 when using the autocorrelation estimator. The results demonstrate that an RCA can, with just 62 receive channels, measure complex 3-D flow fields at a high volume rate.

Pressure Difference Estimation in Non-stenotic Carotid Bifurcation Phantoms Using Vector Flow Imaging.

Local pressure differences estimated using vector flow imaging (VFI) and direct catheterization in seven carotid bifurcation phantoms were compared with simulated pressure fields. VFI correlated strongly with simulated peak pressure differences (r = 0.99, p < 0.00001), with an average overestimation of 12.3 Pa (28.6%). The range between the lowest and highest pressure difference of VFI underestimated simulations by 4.6 Pa (8.06%; r = 0.99, p < 0.0001). The catheter method exhibited no correlation (r = -0.09, p = 0.85). Ten repeated measurements on one phantom revealed a small standard deviation (SD) for VFI (SD = 8.4%, mean estimated SD = 11.5%), but not for the catheter method (SD = 785.6%). An in vivo peak systolic pressure difference of 97.9 Pa (estimated SD = 30%) was measured using VFI in one healthy individual. This study indicates that VFI pressure difference estimation is feasible in phantoms and in vivo and realistic estimates of the SD can be attained from the data.

Fast 3-D Velocity Estimation in 4-D Using a 62 + 62 Row-Column Addressed Array.

This article presents an imaging scheme capable of estimating the full 3-D velocity vector field in a volume using row-column addressed arrays (RCAs) at a high volume rate. A 62 + 62 RCA array is employed with an interleaved synthetic aperture sequence. It contains repeated emissions with rows and columns interleaved with B-mode emissions. The sequence contains 80 emissions in total and can provide continuous volumetric data at a volume rate above 125 Hz. A transverse oscillation cross correlation estimator determines all three velocity components. The approach is investigated using Field II simulations and measurements using a specially built 3-MHz 62 + 62 RCA array connected to the SARUS experimental scanner. Both the B-mode and flow sequences have a penetration depth of 14 cm when measured on a tissue-mimicking phantom (0.5-dB/[ [Formula: see text]] attenuation). Simulations of a parabolic flow in a 12-mm-diameter vessel at a depth of 30 mm, beam-to-flow angle of 90°, and xy-rotation of 45° gave a standard deviation (SD) of (3.3, 3.4, 0.4)% and bias of (-3.3, -3.9, -0.1)%, for ( vx , vy , and vz ). Decreasing the beam-to-flow angle to 60° gave an SD of (8.9, 9.1, 0.8)% and bias of (-7.6, -9.5, -7.2)%, showing a slight increase. Measurements were carried out using a similar setup, and pulsing at 2 kHz yielded comparable results at 90° with an SD of (5.8, 5.5, 1.1)% and bias of (1.4, -6.4, 2.4)%. At 60°, the SD was (5.2, 4.7 1.2)% and bias (-4.6, 6.9, -7.4)%. Results from measurements across all tested settings showed a maximum SD of 6.8% and a maximum bias of 15.8% for a peak velocity of 10 cm/s. A tissue-mimicking phantom with a straight vessel was used to introduce clutter, tissue motion, and pulsating flow. The pulsating velocity magnitude was estimated across ten pulse periods and yielded an SD of 10.9%. The method was capable of estimating transverse flow components precisely but underestimated the flow with small beam-to-flow angles. The sequence provided continuous data in both time and space throughout the volume, allowing for retrospective analysis of the flow. Moreover, B-mode planes can be selected retrospectively anywhere in the volume. This shows that tensor velocity imaging (full 3-D volumetric vector flow imaging) can be estimated in 4-D ( x, y, z, and t ) using only 62 channels in receive, making 4-D volumetric imaging implementable on current scanner hardware.

Flow Complexity Estimation in Dysfunctional Arteriovenous Dialysis Fistulas using Vector Flow Imaging.

Non-invasive assessment is preferred for monitoring arteriovenous dialysis fistulas (AVFs). Vector concentration assesses flow complexity, which may correlate with stenosis severity. We determined whether vector concentration could assess stenosis severity in dysfunctional AVFs. Vector concentration was estimated in four stenotic phantoms at different pulse repetition frequencies. Spectral Doppler peak velocity and vector concentration were measured in 12 patients with dysfunctional AVFs. Additionally, 5 patients underwent digital subtraction angiography (DSA). In phantoms, vector concentration exhibited an inverse relationship with stenosis severity and was less affected by aliasing in severe stenoses. In nine stenoses of 5 patients undergoing DSA, vector concentration correlated strongly with stenosis severity (first stenosis: r = -0.73, p = 0.04; other stenoses; r = -0.69, p = 0.02) and mid-stenotic diameter (first stenosis: r = 0.87, p = 0.006; other stenoses: r = 0.70, p = 0.02) as opposed to peak velocities (p > 0.05). Vector concentration is less affected by aliasing in severe stenoses and correlates with DSA in patients with dysfunctional AVF.

Noninvasive Estimation of Pressure Changes Using 2-D Vector Velocity Ultrasound: An Experimental Study With In Vivo Examples.

A noninvasive method for estimating intravascular pressure changes using 2-D vector velocity is presented. The method was first validated on computational fluid dynamic (CFD) data and with catheter measurements on phantoms. Hereafter, the method was tested in vivo at the carotid bifurcation and at the aortic valve of two healthy volunteers. Ultrasound measurements were performed using the experimental scanner SARUS, in combination with an 8 MHz linear array transducer for experimental scans and a carotid scan, whereas a 3.5-MHz phased array probe was employed for a scan of an aortic valve. Measured 2-D fields of angle-independent vector velocities were obtained using synthetic aperture imaging. Pressure drops from simulated steady flow through six vessel geometries spanning different degrees of diameter narrowing, running from 20%-70%, showed relative biases from 0.35% to 12.06%, depending on the degree of constriction. Phantom measurements were performed on a vessel with the same geometry as the 70% constricted CFD model. The derived pressure drops were compared to pressure drops measured by a clinically used 4F catheter and to a finite-element model. The proposed method showed peak systolic pressure drops of -3 kPa ± 57 Pa, while the catheter and the simulation model showed -5.4 kPa ± 52 Pa and -2.9 kPa, respectively. An in vivo acquisition of 10 s was made at the carotid bifurcation. This produced eight cardiac cycles from where pressure gradients of -227 ± 15 Pa were found. Finally, the aortic valve measurement showed a peak pressure drop of -2.1 kPa over one cardiac cycle. In conclusion, pressure gradients from convective flow changes are detectable using 2-D vector velocity ultrasound.

Accuracy and Precision of a Plane Wave Vector Flow Imaging Method in the Healthy Carotid Artery.

The objective of the study described here was to investigate the accuracy and precision of a plane wave 2-D vector flow imaging (VFI) method in laminar and complex blood flow conditions in the healthy carotid artery. The approach was to study (i) the accuracy for complex flow by comparing the velocity field from a computational fluid dynamics (CFD) simulation to VFI estimates obtained from the scan of an anthropomorphic flow phantom and from an in vivo scan; (ii) the accuracy for laminar unidirectional flow in vivo by comparing peak systolic velocities from VFI with magnetic resonance angiography (MRA); (iii) the precision of VFI estimation in vivo at several evaluation points in the vessels. The carotid artery at the bifurcation was scanned using both fast plane wave ultrasound and MRA in 10 healthy volunteers. The MRA geometry acquired from one of the volunteers was used to fabricate an anthropomorphic flow phantom, which was also scanned using the fast plane wave sequence. The same geometry was used in a CFD simulation to calculate the velocity field. Results indicated that similar flow patterns and vortices were estimated with CFD and VFI in the phantom for the carotid bifurcation. The root-mean-square difference between CFD and VFI was within 0.12 m/s for velocity estimates in the common carotid artery and the internal branch. The root-mean-square difference was 0.17 m/s in the external branch. For the 10 volunteers, the mean difference between VFI and MRA was -0.17 m/s for peak systolic velocities of laminar flow in vivo. The precision in vivo was calculated as the mean standard deviation (SD) of estimates aligned to the heart cycle and was highest in the center of the common carotid artery (SD = 3.6% for velocity magnitudes and 4.5° for angles) and lowest in the external branch and for vortices (SD = 10.2% for velocity magnitudes and 39° for angles). The results indicate that plane wave VFI measures flow precisely and that estimates are in good agreement with a CFD simulation and MRA.

A method to investigate the biomechanical alterations in Perthes' disease by hip joint contact modeling.

Perthes' disease is a destructive hip joint disorder characterized by malformation of the femoral head in young children. While the morphological changes have been widely studied, the biomechanical effects of these changes still need to be further elucidated. The objective of this study was to develop a method to investigate the biomechanical alterations in Perthes' disease by finite element (FE) contact modeling using MRI. The MRI data of a unilateral Perthes' case was obtained to develop the three-dimensional FE model of the hip joint. The stress and contact pressure patterns in the unaffected hip were well distributed. Elevated concentrations of stress and contact pressure were found in the Perthes' hip. The highest femoral cartilage von Mises stress 3.9 MPa and contact pressure 5.3 MPa were found in the Perthes' hip, whereas 2.4 MPa and 4.9 MPa in the healthy hip, respectively. The healthy bone in the femoral head of the Perthes' hip carries additional loads as indicated by the increase of stress levels around the necrotic-healthy bone interface. Identifying the biomechanical changes, such as the location of stress and contact pressure concentrations, is a prerequisite for the preoperative planning to obtain stress relief for the highly stressed areas in the malformed hip. This single-patient study demonstrated that the biomechanical alterations in Perthes' disease can be evaluated individually by patient-specific finite element contact modeling using MRI. A multi-patient study is required to test the strength of the proposed method as a pre-surgery planning tool.